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' Abstract 

We present an efficient method of reconstructing the parameters of the Constrained MSSM 
from assumed future LHC data, applied both on their own right and in combination with the 
cosmological determination of the relic dark matter abundance. Focusing on the ATLAS SU3 
benchmark point, we show that our method can recover the values of its parameters in a prior- 
independent way. We quantify the improvement in constraining the parameters by applying 
lO ' additional information about the relic abundance at the level of WMAP accuracy. We show that 

a further tightening of the cosmological dark matter constraints, as expected from Planck, will 
■ have only a very limited impact. Further external data may be required to break some remaining 

degeneracies. We argue that the method presented here is applicable to a wide class of low- 
energy effective supersymmetric models, as it does not require to deal with purely experimental 
issues, e.g., detector performance, and has the additional advantages of computational efficiency. 
Furthermore, our approach allows one to distinguish the effect of the model's internal structure 
and of the external data on the final parameters constraints. 
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1 Introduction 



If softly broken low-energy supersymmetry (SUSY) provides a correct description of the particle 
physics realm at energy scales around a few hundred GeV and above, then superpartners are likely 
to be discovered at the LHC. One of the main goals of ATLAS and CMS experiments will be to 
identify those particles by determining their masses and other properties. 

The actual outcome will depend not only on the LHC machine and detector performance but 
obviously also on the mass scales of the superpartners themselves. A whole plethora of different 
possibilities can be listed here, ranging from one extreme where all of the superpartners may come 
out to be too heavy to for the LHC reach, to another where all, or most, of them will be discov- 
ered. Unfortunately, basically the whole spectrum of options remains open even in perhaps the 
most economical SUSY framework, the Constrained Minimal Supersymmetric Model (Constrained 
MSSM, or CMSSM) [lj which includes the minimal supergravity (mSUGRA) model [2], as shown 
by a number of recent global fits of the CMSSM based on Bayesian statistics [HI 0] and on a % 2 
approach 16] . While the latter show a stronger preference for a fairly low SUSY mass scale, in the 
range of a few hundred GeV, the former point to a more cautious picture, where a much wider mass 
range remains allowed. This discrepancy is caused by the fact that, with the data that is currently 
available, even the CMSSM still remains to some extent underconstrained, and the specifics of the 
statistical and data analysis treatment can lead to fairly different results. It is therefore clear that 
selecting, or at least limiting, SUSY models by using LHC measurements is certainly going to be a 
very challenging task as there exist large degeneracies among the MSSM parameters that can lead 
to indistinguishable LHC signatures (see, e.g., Ref. [7]). 

In preparation for dealing with real data, a number of approaches to particle mass reconstruction 
have been developed based on extracting kinematic information from one or more decay chains of 
superpartners, typically requiring two or more visible particles in the final state [8]. These and other 
techniques have been used by LHC experimental groups which have performed a large number of 
detailed studies in a few reference, or "benchmark", points, often selected in such a way as to 
typically allow several of the superpartners to be seen at the LHC. 

In a recent extensive ATLAS Report [9], in the framework of the CMSSM/mSUGRA a so-called 
ATLAS SU3 benchmark point (which is specified below) was examined with Markov Chain Monte 
Carlo (MCMC) scans, with the aim of evaluating the expected accuracy of reconstructing CMSSM 
parameters: a common gaugino mass parameter m^j a common scalar mass parameter mo, a 
common trilinear term Aq, all evaluated at the unification scale Mgut — 2 x 10 16 GeV, plus a 
Higgs vacuum expectation values tan/3. Assuming an integrated luminosity of lfb -1 , a dilepton 
and lepton+jets edge analysis of the decay chain qi —>■ Xq,{~^ I l^)q — > Xi^ + ^Q an d the high-p^ 
and large missing energy analysis of the decay chain qa — > XiQ were performed, where qi (qn) 
denotes the first or second generation left (right) squark, Xi 2 the first and second neutralino and I 
an intermediate slepton. It was concluded that m 1( / 2 and mo could be reconstructed with adequate 
accuracy, while prospects for tan/3 looked somewhat poorer, and even more so for Aq (see Ref. [9], 
pp. 1617, ff). 

In this paper we perform an independent analysis of the ATLAS SU3 point using the publicly 
available information about the expected ATLAS capabilities to measure the SU3 mass spectrum. 
We first demonstrate that a simple modelling of the mass spectrum constraints in an effective 
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likelihood is sufficient to reproduce with reasonable accuracy the results of the full ATLAS analysis, 
while being much more economical in terms of computational requirements. We then build on the 
ATLAS analysis by examining in detail the impact of priors on the induced constraints on the 
CMSSM parameters, demonstrating that, with expected LHC data, one will be able to conclusively 
resolve the prior dependence which currently affects analyses employing presently available data. We 
also compare with the limits that can be obtained using a maximum likelihood analysis, and we show 
that the choice of statistics (Bayesian or maximum likelihood) no longer matter once one combines 
ATLAS data with cosmological relic abundance determinations. We clarify what role played by 
assuming a specific theoretical model (here the CMSSM) in complementing the information coming 
from the ATLAS measurements with model-specific theoretical correlations between masses of the 
observables. 

In the second step, we go beyond the ATLAS analysis by applying additional information about 
the cosmological relic density Q x h 2 of the lighest neutralino x? (below often denoted by x f° r 
simplicity), assumed to be dark matter (DM) in the Universe. The dark matter abundance clearly 
provides additional information about the model at hand, and in this analysis we aim at obtaining a 
quantitative measure for the extra constraining power that it provides, on top of that expected from 
the ATLAS data. Here we first impose WMAP uncertainties on Q x h 2 and demonstrate a significant 
improvement in the determination of the CMSSM parameters, especially thq. Next we investigate 
the impact of further reducing the observational errors on Q x h 2 to an accuracy as expected from 
Planck, and show that this will lead to only a further modest improvement. We also comment on 
the impact of some other commonly used constraints, in particular from b — > 57 and {g — 2)^, which 
we, however, do not apply here. Finally, we examine the impact of the different constraints from 
ATLAS and cosmology on the uncertainties of mass measurements of several superpartners and on 
predictions for the scattering cross section relevant for direct detection of dark matter experiments. 

The paper is organized as follows. In Sec. [2] we present a setup of our analysis for the ATLAS 
SU3 point in the CMSSM, including details of our scans of the CMSSM parameters for two different 
priors. In Sec. [3] we present our numerical results for the posterior probability density functions 
(pdfs) demonstrating their prior-independence, and in SecH]we compare them with the alternative 
measure of the profile likelihood. We summarize our findings in Sec. 

2 Setup and benchmark point 
2.1 The ATLAS SU3 benchmark point 

We examine the ATLAS SU3 benchmark point for which input values of CMSSM parameters are 
given on the left side of Table [TJ Since in the ATLAS analysis errors of relevant SM parameters 
("nuisance parameters") were not included, we assume that the benchmark values for the nuisance 
parameters are set at their central values as given on the right side of Table [U In the reconstruc- 
tion done below, we then allow the nuisance parameters to vary and we constrain them using the 
likelihood given below. 

Since LHC data is rather unlikely to differentiate among the flavors of the squarks of the first 
two generation, in what follows we denote them all by a common symbol q, and by m,g their average 
mass, similarly as in Ref. [9]. On the other hand, I will denote the lightest slepton and my- its mass. 
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CMSSM parameter ATLAS SU3 benchmark value 


SM parameter Input value 


m 1/2 300 GeV 
m 100 GeV 
tan/3 6.0 
Aq -300 GeV 
sgn(^) + 


M t 172.6 GeV 
mbimb) 1 ™ 4.20 GeV 
asiMz)^ 0.1176 
l/ae^Mz) 7 ^ 127.955 



Table 1: Left side: input CMSSM parameters values for the ATLAS SU3 benchmark point. Right 
side: input values of relevant SM parameters used in the numerical analysis. 



superpartner 


mass 


superpartner 


mass 


superpartner 


mass 


xi(= x) 


117.9 GeV 




230.8 GeV 




666.2 GeV 


X°2 


223.4 GeV 


en, Jj-R 


157.5 GeV 


da 


639.0 GeV 


x° 3 


463.8 GeV 




217.5 GeV 


UL 


660.3 GeV 


xl 


479.9 GeV 


n 


152.2 GeV 




644.3 GeV 


xt 


224.4 GeV 




232.4 GeV 


h 


599.0 GeV 


xt 


476.4 GeV 


V T 


216.9 GeV 


b 2 


636.6 GeV 


9 


717.5 GeV 






h 


446.9 GeV 










t 2 


670.9 GeV 



Table 2: Superpartner mass spectrum for the ATLAS SU3 point. 



In the case of the ATLAS SU3 point its role is played by t\. 

The resulting mass spectrum, as computed using the SoftSusy code [10] in the 1-loop approx- 
imation is given in Table By comparing with the mass spectrum for the ATLAS SU3 point 
given in Ref. [9] (cf. Table 2 on page 1516), we can see some differences, especially a systematic 
shift in squark masses by a few tens of GeV, which may be due to using different numerical codes, 
approximations (although in both cases 1 loop expressions are applied) as well as different inputs in 
SM parameter values. While this will contribute to some differencies we will find with the ATLAS 
results, at the end these discrepancies are of secondary importance, as we discuss below. 

2.2 The likelihood function 

The study performed by the ATLAS Collaboration on the SU3 point reports the expected accuracy 
in the reconstruction of some of the masses and mass differences in the benchmark SUSY spectrum 
given in Table [2j Dilepton edges will constrain m^o, m x o, m-q and my with fairly poor accuracy, while 
providing much tighter limits on the mass differences between the three latter quantities and the 
lightest neutralino, since these follow more directly from endpoint measurements. The observable 
quantities to be constrained by ATLAS are given by the set 

8 = {m x o,m x o — m x o,mj— m x »,mj — m x o}. (1) 

We further assume that the maximum likelihood (ML) value of 8 obtained by ATLAS, 6ml, cor- 
responds to the value of the true benchmark point, 6* = {117.9, 105.5, 34.3, 534.5} GeV, where the 
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m^o m^o — m v o mr — m^o m?r — m^o 

Al A3 Al t Al y Al 


m v o 

m^o — m^o 

A2 Al 

m r -m x o 


3.72 x 10 3 53.40 1.92 x 10 3 10.75 x 10 2 
3.6 29.0 -1.3 
1.12 x 10 3 4.65 
14.1 



Table 3: ATLAS covariance matrix employed in the analysis. 



numerical value is obtained from Table [2j In other words, we neglect realization noise, an assump- 
tion which is justified by the fact that {6ml) = 9*, where (•) denotes an average over realizations. 
The likelihood function from ATLAS is then modeled as a Gaussian centered around the true value 
of the observable quantities, 

- 2 In £atlas = Xatlas = (0 - 0ml)' C- 1 (9 - 9ml), (2) 

where the covariance matrix C is given in Table [3] (and we have dropped an irrelevant normalization 
constant). It represents the full covariance between the masses and the mass differences. The 
covariance matrix includes statistical errors only; systematic errors are negligible. 

The form of the ATLAS likelihood function given in Eq. ([2]) is a simple Gaussian approximation 
to the actual likelihood function that one would obtain from a full analysis of simulated ATLAS data. 
The latter is, however, not available outside the Collaboration, and therefore our approximation 
represents the best that can be reasonably done given the information that is expected to be publicly 
available. There are two reasons why it might be interesting to consider an approximate ATLAS 
likelihood function at the level of the SUSY mass spectrum. Firstly, it is not unreasonable that the 
simple approximation adopted here will give a fairly accurate representation of ATLAS capabilities 
(see below for further comments on how this compares with the full analysis carried out by the 
ATLAS Collaboration), and therefore provide a useful shortcut to quantitatively implementing the 
constraints that would result from a full analysis. Secondly, given the model-independent constraints 
on the low-energy SUSY mass spectrum, one can use them to constrain the high-energy parameters 
of any SUSY model of the MSSM class, in this case the CMSSM. As we shall show below, the 
constraining power included in the low-energy likelihood is then supplemented by the theoretical 
structure of the model itself. 

We will now examine prospects for reconstructing the input values of the ATLAS SU3 point. 
We consider the following data combinations: 

• ATLAS data only: including only the likelihood function given by Eq. ([2]); 

• ATLAS data (as above) + WMAP-level uncertainty on dark matter abundance; 

• ATLAS data (as above) + Planck-level uncertainty on dark matter abundance. 

When including information about the cosmological dark matter abundance, we have to address 
the fact that the value of the neutralino relic abundance (computed using the code Micromegas [TT] ) 
for the ATLAS SU3 benchmark point, is Q. x h 2 = 0.2332, which is a factor of some 2.5 above the 
WMAP range of 0.1099 ±0.0062 for the cosmological dark matter |12j . However, assuming standard 
Big Bang cosmology and that the CMSSM is correct, we would expect that for the actual measured 
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value of the CMSSM parameters, for which ATLAS and/or CMS measurements are made, VL x h? 
would lie in the WMAP range. It is therefore not unreasonable to examine the impact of the 
extra piece of information, both in terms of its value and uncertainty, which is provided by the 
cosmological relic density of dark matter. Notice that we do not rely here on a potential ability to 
reconstruct Q x h 2 from LHC data alone [13] (which, on the other hand, could provide an important 
cross-check of our cosmological input) but use it as external constraint. Therefore, since the ATLAS 
analysis has been carried out for the ATLAS SU3 point, with the specific values of the CMSSM 
parameters, we also adopt to keep this central value for the relic abundance. On the other hand, 
we don't expect any major change in the numerical results presented below if Q x h 2 were actually 
close to the WMAP range. 

In our analysis we implement the cosmological dark matter determination as follows. Firstly, 
we assume that the true value of the relic abundance corresponds to the value computed for 
the benchmark point, and that WMAP-level constraints correspond to a Gaussian likelihood cen- 
tered around that value with standard deviation given by the current WMAP uncertainty, namely 
cwmap = 6.2 x 1(T 3 : 

- 21n£ W MAP = Xwmap = - — — 2 — ■ ( 3 ) 

a WMAP 

The Planck satellite is expected to improve the WMAP accuracy on the relic abundance by a factor 
~ 10 |14j . When including Planck-level constraints, we therefore adopt the same likelihood as above 
but with a smaller standard deviation, <7pi a nck = 6.2 x 10 . Since for the "bulk region" to which the 
ATLAS SU3 point belongs the theoretical error in the relic abundance is estimated to be tiny [15] 
we neglect it here. In other cases it can be much larger, primarily due to the larger uncertainties 
in computing mass spectra; for example in the focus point region it would likely dominate and this 
would cloud the potential impact of the cosmological data. 

In order to facilitate a comparison with the ATLAS study [9j, in this analysis we do not apply 
any other constraints, e.g., from LEP, rare processes in heavy quark physics (6 — ► sy, etc) or the 
anomalous magnetic moment of the muon (g — 2)^, which are routinely used in global analyses 
of the CMSSM and other popular SUSY models. We have also checked that fixing the nuisance 
parameters or marginalising over them has a negligible impact on the results. Therefore we only 
present results with the nuisance SM parameters fully marginalised, even though in the ATLAS 
analysis the SM parameters were fixed at their central values. 



2.3 Scanning the CMSSM parameters 

With the aim of reconstructing the true values of the defining parameters for the ATLAS SU3 
point, we explore the CMSSM parameter space with the help of the Nested Sampling (NS) scanning 
technique, as implemented in the MultiNest algorithm |16j . 
We consider two different priors: 

• flat prior: flat in mi/2, mo, Aq, tan/3, with the ranges: 50 GeV < mj^jmo < 500 GeV, 
2 < tan (3 < 62 and -4 TeV < Aq < 4 TeV. 

• log prior: flat in log mi/ 2 , log tuq, Aq, tan 0, with the same ranges: log(50) < log mo, log mi/ 2 < 
log(500) (in GeV), and as above for Aq and tan/3. 
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Notice that we have employed here narrower ranges of miw and mo than the values of up to 
a few TeV used in our previous analyses [U Q2]. However, we have checked that enlarging the 
prior range to much larger values of m\/2,mo (up to 4 TeV) has no impact on our reconstructed 
parameter values, as our algorithm correctly recovers the true parameter values even in the case of 
a much larger prior range. Finally, for the SM nuisance parameters we assume the same ranges as 
in our previous papers [HUB]; in any case as we have mentioned above, the details of the treatment 
of nuisance parameters has basically no impact on the results presented here. 

3 Results 

In this Section we present our numerical results from scans performed using a modified version of 
the publicly available SuperBayes package [19]. 



Roszkowski, Ruiz de Austri & Trotta (2009) Roszkowski, Ruiz de Austri & Trotta (2009) 




rri (GeV) A n (TeV) 



Figure 1: 2D posterior pdf for the case of applying ATLAS mass spectrum data alone, for some 
CMSSM parameter combinations and for the log prior choice. Compare with Fig. 12, p. 1638 of [9]. 



3.1 Using ATLAS data only 

We begin by considering the constraining power on the CMSSM parameters of ATLAS mass spec- 
trum data alone. In Figs. Q] and [2] we present 2-dimensional (2D) Bayesian posterior pdfs assuming 
the log and the flat prior, respectively, while the corresponding ID pdfs for the log prior case are 
shown in Fig. [3J By examining the (m^2,mo) plane for the log prior in Fig. [TJ we can see that 
both the 68% (inner regions) and the 95% (outer regions) total probability regions are well centered 



6 




Figure 2: The same as in Fig. Q] but for the flat prior case. Notice that the two choices of priors 
lead to statistically indistinguishable results. 

around the true value. Our algorithm recovers the best-fit point within the limits of numerical noise 
for all CMSSM parameters, except for Aq, where the reconstructed best-fit ends up in a wrong re- 
gion of the parameter space due to a fundamental degeneracy discussed below. The posterior mean 
is also reasonably close to the true value (within la for all parameters, except for Aq), although it 
is slightly skewed due to the asymmetric nature of the contours, which exhibit heavier tails than 
Gaussian (see also Fig. [3] below). On the other hand, tan/3 is somewhat less well reconstructed, 
yielding only an upper limit. 

In contrast, Aq is rather poorly constrained in this case, and actually shows a sign ambiguity. 
This is because it enters the analysis in a rather indirect way, mostly via the off-diagonal terms 
X T = A T — [i tan (5 in the stau mass matrix, where A T is the value of Aq at the EW scale evaluated 
with its RGE and /i is computed from the usual requirement of correct electroweak symmetry 
breaking. A closer examination reveals that, for Aq ~ 1 TeV (in between the two la regions in 
the right panel of Fig. []]), X T is minimized and the mass difference between T2 and t\ (which plays 
the role of the lightest slepton in the decay chain) goes to zero. Since in the ATLAS analysis only 
r\ was^ considered, such cases are not allowed. Our study thus reveals that in studying the decay 
X2 —> I F the exchange of both T2 and t\ should be considered, as for some values of Aq 

their masses, and therefore also relative contributions, may be comparable. 

In the case of the flat prior (Fig. [2]) the emerging picture remains essentially identical, thus 
confirming that the prior choice becomes irrelevant once the constraining power of the data is 
sufficiently strong, as expected. Many of the features seen in Fig. [1] are displayed more clearly in 
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Fig. [3] where the corresponding ID pdfs are presented for the log prior case only; the flat prior 
produces very similar results and is therefore not shown. We give the 68% and 95% intervals of our 
reconstructed CMSSM parameters in Table 01 A comparison with the profile likelihood is carried 
out further below. 



Applying ATLAS data only 


Parameter 


True value 


Best fit 


68% (95%) range 


m i/2 (GeV) 
m (GeV) 
tan (3 
A (GeV) 


300 
100 
6.0 
-300 


300.4 
104.3 
3.8 

1749.7 


[288.2,308.4], ([278.2,316.3]) 
[98.7,173.6], ([89.8,235.2]) 
< 13.8 (< 27.4) (1 tail) 
[-568.3, 1701.8], ([-995.1, 2311.6]) 


m x o (GeV) 

n x k 2 

logo£ J (pb) 


117.9 

0.2332 

-8.92 


116.9 

0.2330 

-8.87 


[113.7,120.8], ([110.6,123.7]) 
[0.2264,0.2870], ([0.2096,0.3450]) 
[-9.14,-8.42], ([-9.45,-8.04]) 


Applying ATLAS+WMAP-like data 


Parameter 


True value 


Best fit 


68% (95%) range 


m i/2 (GeV) 
m (GeV) 
tan (3 
A (GeV) 


300 
100 
6.0 
-300 


302.3 

98.3 

5.5 

-228.2 


[293.2,310.7], ([285.5,317.5]) 

[95.9,112.2], ([90.9,151.6]) 

< 7.3 (< 16.3) (1 tail) 

[-498.1, 1437.6], ([-887.7, 2199.1]) 


m x o (GeV) 

n x h 2 

log (pb) 


117.9 

0.2332 

-8.92 


118.6 

0.2333 

-8.85 


[115.1,121.3], ([112.2,123.8]) 
[0.2281,0.2397], ([0.2225,0.2454]) 
[-9.07,-8.51], ([-9.36,-8.03]) 


Applying ATLAS +Planck-like data 


m i/2 (GeV) 
m (GeV) 
tan/3 
A (GeV) 


300 
100 
6.0 
-300 


300.5 

99.4 

6.1 

-257.4 


[295.7,311.1], ([289.0,317.6]) 

[95.3,106.1], ([92.0,115.6]) 

< 4.3 (< 11.3) (1 tail) 

[-397.5, 1378.7], ([-700.1, 2045.5]) 


m x o (GeV) 

n x h 2 

log of 7 (pb) 


117.9 

0.2332 

-8.92 


118.0 

0.2332 

-8.88 


[115.9,121.3], ([113.3,123.8]) 
[0.2327,0.2338], ([0.2322,0.2345]) 
[-8.99,-8.56], ([-9.20,-8.31]) 



Table 4: Reconstructed values and errors for the input CMSSM parameters and for some key 
observables. We also give the best fit from our scan. The 68% and 95% ranges are computed from 
the posterior pdf as shortest intervals around the mean. For definiteness, we have employed the log 
prior scan but the results from the flat prior case are essentially identical. 

When considering posterior constraints on the SUSY mass spectrum, it is apparent that some 
of the constraints are much stronger than the likelihood function alone (which actually applies 
to a more general case of the MSSM) would seem to imply. For instance, the la error on m x o 
from the ATLAS likelihood is 60 GeV (cf. Table El where the likelihood la range is obtained as 
the square root of the diagonal elements). However, the reconstructed neutralino mass within 
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Figure 3: ID posteriod pdf for the case of applying only ATLAS mass spectrum data, summarized 
in ID projections. We show only the log prior case, for the flat prior case is essentially identical. It 
is clear that ATLAS data alone is not sufficient to reconstruct all of the CMSSM parameters. In 
particular, while mi/2 is well measured, Aq and tan/3 remain largely undetermined. 
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Figure 4: Illustration of the extra constraints coming from the assumption of the CMSSM as the 
theoretical framework. The parameter space accessible within the CMSSM is given by the green 
region (dots represents uniformly drawn samples), the red ellipses are the ATLAS likelihood (as 
given by the covariance matrix of Table E]) while the blue contours are the posterior constraints. In 
the context of the CMSSM this allows to derive much tighter constraints on m o than it would be 
possible based on the likelihood alone. 



the CMSSM shows a much smaller error, of order ~ 4 GeV, cf. Table [H The reason for this is 
that the information supplied by the likelihood is supplemented by the internal structure of the 
CMSSM parameter space, within which the masses of many of the sparticles are highly correlated. 
This is demonstrated in Fig. HJ where one can see that the correlation between masses in the 
spectrum within the CMSSM is nearly orthogonal to the constraints provided by ATLAS for the 
mass spectrum observables plotted in the Figures. One can think of this correlation as an additional 
a priori piece of information contained in the model. In other words, given the theoretical structure 
of the CMSSM, certain mass combinations in the spectrum (which are otherwise allowed by the 
projected ATLAS constraints, red ellipses in Fig. H]) are simply not allowed by the structure of the 
model. Therefore the final constraints on the spectrum are much tighter than the likelihood alone 
would imply. Supplementing the mass spectrum constraints with a model-specific implementation, 
as done here, has the additional advantage of displaying which part of the constraining power comes 
from the experimental data and which one from the theoretical properties of the model. 

It is interesting to examine how well ATLAS data alone can determine the relic abundance for 
the ATLAS SU3 point, in some analogy with what, for example, has been done for some other 
benchmark points in Refs. [13} I20j . This is shown in Fig. [5] where we find that, from the assumed 
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ATLAS data alone one would obtain Q x h 2 = 0.253 ± 0.034, hence with a relative accuracy of 
~ 13%. For this specific point, this would imply that ATLAS data would determine the neutralino 
dark matter abundance at about la. Since the neutralino dark matter abundance for the ATLAS 
SU3 benchmark point is some 2.5 times larger than the value currently preferred by cosmological 
observations, if we assume that the accuracy for the benchmark point is representative for the 
accuracy that ATLAS will actually find around a point with the correct cosmological relic abundance 
of about 0.11, our estimate is that ATLAS data alone would be able to determine the DM relic 
abundance at the level of ~ 3a. Finally, in our present analysis we have ignored any theoretical 
error in the DM abundance prediction. While for the ATLAS SU3 point, which falls into the bulk 
region, such an error is likely to be very small, in general it should be folded in when producing the 
posterior shown in Fig. [5j 

Roszkowski, Ruiz de Austri S Trotta (2009) 
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Figure 5: ID posteriod pdf for the relic DM abundance f^/i 2 of the neutralino, obtained from 
ATLAS data alone. We show only the log prior case, for the flat prior case is basically the same. 

3.2 Comparison with the ATLAS Collaboration results 

It is instructive to compare Figs. [1] and [2] with Fig. 12 in Section 9.3 of the ATLAS Report [9] 
where 2D marginal Bayesian posteriors are presented following the ATLAS MCMC analysisEI The 

4 Although it is labeled "likelihood maps", the quantity plotted in Fig. 12 in of Ref. [9] is actually a marginal 
Bayesian posterior (Peter Wienemann - Private Communication), analogous to the one plotted in our Figs. [T] and [2] 




68%, 95% contours 
® best fit 
mean 
^ true value 
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overall shape of the high-probability (m^,^) an d (tan/3, Ao) regions is qualitatively similar al- 
though quantitatively we find somewhat less stringent bounds. In particular, we can see the largest 
difference in the case of Aq where the highest probability region found in [9] lies on the boundary 
of the correct region found in our analysis, while the other, multi-TeV region, is in Ref. [U] absent 
altogether. There is also some difference in tuq which in our case is not as well constrained as in 
Ref. 0. 

It is however difficult to carry out a closer comparison, since not many details are given regarding 
the setup used in the ATLAS fitting analysis, in particular, about their treatment of SM nuisance 
parameters. Also, the ATLAS fitting analysis was performed directly from end-point measurements 
while we used a Gaussian approximation to the likelihood for masses and mass differencies alone, 
thus inevitably loosing a certain amount of information contained in the full analysis. It is, however, 
certainly encouraging that our "shortcut" method of reconstructing CMSSM parameters using a 
relatively crude approximation to the full ATLAS analysis was able to recover quite compatible 
regions of SUSY parameters around their true values. The only exception is Aq, as explained 
above. As we show below, adding cosmological relic abundance constraints does help in further 
tightening some of the constraints. We conclude that, despite those differences, overall we find a 
reasonably good agreement with the ATLAS analysis. This suggests that not too much information 
is lost by carrying out the analysis employing an effective likelihood at mass spectrum level. 

The advantage of our procedure is that it allows one to easily change the model-specific assump- 
tions: if one replaces the CMSSM by another SUSY model that one is interested in, the analysis can 
be carried out without the need of going through the details of detector performance and obtaining 
the ATLAS likelihood numerically via Monte Carlo, thereby strongly reducing the computational 
requirements. In fact, our analysis requires about 24 hours on eight 3GHz processors, and it is 
therefore relatively computationally undemanding. Furthermore, it would be easy to adapt our 
method to employ a more complete likelihood function on the mass spectrum should this become 
available as part of the data products released by the LHC Collaborations. This would allow theo- 
retical studies of the constraints implied on different SUSY models without the need to reproduce 
the full detector-specific signal reconstruction. 

3.3 Impact of including the DM relic abundance 

We now add to our likelihood function a constraint on the relic abundance uncertainties, as discussed 
in Sec.EJ In Fig.[6]we show the effect of imposing the ATLAS and WMAP data (ATLAS+WMAP) 
in the upper row, and an analogous case for the ATLAS+Planck case in the lower row. We plot 
the posterior for the log prior case; the flat prior case is basically identical. It is clear that, adding 
WMAP-like constraints improves the reconstructing power in determining the CMSSM parameters 
very considerably in the case of mo (and to some extent also tan/3), while the impact on the other 
two CMSSM parameters is fairly limited. This can be traced back to the fact that, in the bulk 
region, Q x h 2 is determined primarily by the mass of the lightest slepton, via a i-channel exchange. 
Tightening the allowed range of S^ x /i 2 selects a more peaked range of m^-and thus also mo on which 
it mostly depends. On the other hand, m\i2, which primarily determines m x , can be adequately 
constrained already by using only ATLAS data. As regards Aq, the bi-modality still remains as 
it is caused by the internal structure of the CMSSM. On the other hand, a further improvement 
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Figure 6: Impact of adding to the ATLAS data cosmological dark matter abundance determination 
with WMAP-like (upper row) and Planck-like (lower row) errors on (mi/2, mo) (left panels) and 
(tan/3, Aq) (right panels). Only the log prior case is presented; the flat one produces very similar 
results. 



of the error on Q x h? to the level expected from Planck does not seem to improve the situation 
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much beyond the ATLAS+WMAP case. In this context we again emphasize that, at this level of 
accuracy, it will be essential to achieve a similar, or better, level of theoretical errors, which may 
be challenging even for the bulk region. 

In Fig. [7] we show the constraints on the masses of several superpartners obtainable with the 
three sets of data considered in this paper. We can see that in the case of the gauginos (x?2' Xi 
and g), whose masses are determined primarily by m^ 2 (which is well reconstructed), the errors 
are rather small, while for higgsino-like states (X34 an d X2 ) the errors are large because of a poor 
determination of the jj, parameter. For the states whose mass strongly depends on tjiq (spin-zero 
superpartners) the errors again reflect that of the common scalar mass, whose reconstruction, while 
reasonable, is not as good as for m^. 

Roszkowski, Ruiz de Austri Sc Trotta (2009) 

- ATLAS alone CMSSM, m> - 

800 _ ATLAS+WMAP ATLAS SU3 point _ 

ATLAS+Planck 




Figure 7: Reconstruction of the SUSY mass spectrum using projected ATLAS data only (red), 
adding WMAP-like constraints on the dark matter relic abundance (blue) and adding Planck-like 
dark matter constraints (green). 

The impact of further imposing other often used constraints from b — > s~f and {g — 2)„ is in the 
present case rather limited. This is because the total error in the first quantity is still substantial 
while SUSY contributions to (g — 2)^ for the low mass region are large enough to contribute to the 
apparent discrepancy between the experimental data and the SM value. 

Finally, we investigate how well one can predict the spin-independent cross section a^ 1 of dark 
matter neutralino scattering off a proton tested in direct detection experiments. As can be seen 
from Fig. [8j at 68% the value of 1 will remain uncertain to within about one order of magnitude, 
while the neutralino mass will be very well constrained by LHC data as a reflection of the bounds on 
m-x/2- This is because in the case studied here is too a large extent determined by a i-channel 
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heavy scalar Higgs exchange, where, in addition to m x , the main two parameters are tan/3 and the 
Higgs mass which shows a considerable spread of values, mostly due to the larger uncertainty in 
m,Q. Adding information about the DM relic abundance therefore improves the situation only in a 
fairly limited way. 



-6.5 
-7 
-7.5 
_X -8 

S-B.5 
o 



Roszkowski, Ruiz de Austri & Trotta (2009) 



_Q 
Q. 



-9: 
-9.5: 



68%, 95% contours 
Black: ATLAS only 
Blue: ATLAS+WMAP 
^ true value 



Posterior pdf 
Log priors. 
CMSSM, li>0 
ATLAS SU3 point- 




110 115 120 125 
m x (GeV) 



-6.5 
-7 
-7.5 
_X -8 

C/3 Q_ 

* -8-5 
o 



Roszkowski, Ruiz de Austri & Trotta (2009) 



_Q 
Q. 



-9 
-9.5 



68%, 95% contours 
Black: ATLAS only 
Blue: ATLAS+Planck 
^ true value 



Posterior pdf 
Log priors . 
CMSSM, n>0 
ATLAS SU3 point- 




110 115 120 125 
m x (GeV) 



Figure 8: Impact of adding to the ATLAS data cosmological dark matter abundance determination 
with WMAP-like (left panel) and Planck-like (right panel) errors on the predictions for the spin- 
independent cross section of dark matter neutralino scattering off a proton, relevant for direct 
detection experiments. 



4 Comparison with profile likelihood 

In order to examine the robustness of the results obtained with the Bayesian posterior pdf, in 
this Section we compare them with what would be obtained by using a more traditional x 2 -based 
analysis. We thus define the profile likelihood for, e.g., the CMSSM parameter mi, where mi, . . . , ms 
comprise the 4 CMSSM parameters and the 4 SM nuisance parameters of Table [TJ as 

£(mi) = max C(d\m), (4) 

rn2,...,rns 

where in our case £{d\m) is the full likelihood function. Thus in the profile likelihood one maximises 
the value of the likelihood along the other CMSSM and SM parameters, rather than integrating it 
out as in the marginal posterior. From the profile likelihood, confidence intervals are then obtained 
using the usual likelihood-ratio criterion. In the context of MCMC scans of the parameter space, the 
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Figure 9: Comparison between the profile likelihood (filled) and posterior (empty contours) using 
ATLAS mass spectrum data only (top row) and adding Planck (bottom row). With this combination 
of data, the choice of statistics (Bayesian posterior or profile likelihood) becomes almost irrelevant, 
giving intervals that match at the 10% level. 



profile likelihood can be evaluated by simply finding the maximum likelihood value within a given 
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bin. This has been studied before in the context of the CMSSM in Refs. |17[ 118] . Its interest lies in 
the fact that it is a prior independent measure. Since we have found above little prior dependence 
of the posterior pdf, this suggests that the posterior pdf is dominated by the likelihood. Therefore 
we generically expect that the profile likelihood will give similar statistical results as the posterior 
studied above. 

This is indeed confirmed in the top row of Fig. where we present the 68% (inner contours) and 
the 95% confidence regions (outer contours) of the profile likelihood for the ATLAS-only data case 
and we compare them with the analogous regions derived from the posterior pdf presented earlier 
in Fig. [TJ We show only the log prior case as the profile likelihood is prior independent (which 
we have verified numerically). We can see that at the 95% CL from the profile likelihood is quite 
similar to the corresponding 95% region derived from the Bayesian posterior for all parameters, 
except for mo, for which the profile likelihood yields looser constraint. It is worth noticing that 
the posterior pdf yields a somewhat better reconstruction of m,Q and tan f3 and a similar one for 
the other CMSSM parameters. However, the relative merits of the reconstructed confidence regions 
from the posterior or from the profile likelihood cannot be assessed here. It is in general a difficult 
task to decide which statistics yields the "best" results (however one chooses to define this). A 
possible way forward would be to carry out a coverage study of the quoted confidence intervals, 
which is beyond the scope of this paper. 

On the other hand, what is encouraging is that, when the data becomes sufficiently constraining, 
both statistical quantities produce essentially equivalent confidence intervals. This is presented in 
Fig. [9] for the ATLAS+Planck case, which should be compared with the bottom row of Fig. [HI 

5 Summary and conclusions 

In this paper we have examined prospects for reconstructing supersymmetric parameters from as- 
sumed future data that one can reasonably expect to become available. To this end we focused on 
the Constrained MSSM and on the benchmark point ATLAS SU3. 

By following the ATLAS assumptions as closely as possible without having access to the full 
simulated likelihood function, we arrived at generally rather similar results for the reconstruction of 
the CMSSM parameters, with the exception of Aq, for which our projected limit appears somewhat 
weaker. We highlighted the computational advantage of our method which employs an effective 
likelihood at the mass spectrum level, which allows to shortcut the computationally expensive 
simulation of the whole experimental setup. We also demonstrated that, once LHC data become 
available, previously observed prior dependence of the results disappears. We showed that the 
conclusions depend only mildly on which statistical quantity one chooses to adopt, i.e., Bayesian 
posterior or profile likelihood, in marked contrast with the present-day situation. 

We then extended the analysis by adding to the likelihood function information about the 
neutralino dark matter relic abundance by imposing WMAP-like and Planck-like constraints. This 
improved the ability to reconstruct the value of especially mo and tan/3, much less so for 
(compared to the ATLAS data only case), while the bi-modality in the determination of Aq could 
not be removed. 

While the ATLAS SU3 point (and maybe also the CMSSM in the first place) may be unlikely 
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to be realized in Nature, the method presented here appears to be powerful and robust enough 
to adequately reconstruct supersymmetric parameters from summary statistics of LHC measure- 
ments. The additional advantages presented here are the ability to easily investigate several different 
theoretical scenarios with relatively little computational effort, and the capability to produce pre- 
dictions for derived observable quantities, such as for example the cosmological relic abundance and 
direct detection cross sections. The inclusion of observational constraints from such probes has also 
been demonstrated to be easily implemented. Finally, the favourable scalability of our MultiNest 
scanning algorithm with the dimensionality of the parameter space means that this method is in 
principle ready to investigate theories with several tens of free parameters, thereby opening the way 
to massive inference in supersymmetry phenomenology. As such we believe that our method will 
be a useful tool to face the real data that is expected to soon start arriving from the LHC, even if 
differs significantly from the case considered here. 
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